Non-local Andreev reflection in superconducting quantum dots 
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With the aid of the Keldysh technique we develop a microscopic theory of non-local electron 
transport in three-terminal NSN structures consisting of a chaotic superconducting quantum dot 
attached to one superconducting and two normal electrodes. Our theory fully accounts for non- 
equilibrium effects and disorder in a superconducting terminal. We go beyond perturbation theory 
in tunneling and derive a general expression for the system conductance matrix which remains valid 
in both weak and strong tunneling limits. We demonstrate that the proximity effect yields a decrease 
of crossed Andreev reflection (CAR). Beyond weak tunneling limit the contribution of CAR to the 
non-local conductance does not cancel that of direct electron transfer between two normal terminals. 
We argue that temperature dependence of the non-local resistance of NSN devices is determined by 
the two competing processes - Andreev reflection and charge imbalance - and it has a pronounced 
peak occurring at the crossover between these two processes. This behavior is in a good agreement 
with recent experimental observations. 
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I. INTRODUCTION 

Non-local (crossed) Andreev reflectionii^ is the pro- 
cess which occurs in multi-terminal hybrid normal metal- 
superconductor-normal metal (NSN) proximity struc- 
tures and involves two subgap electrons entering a su- 
perconductor from two different normal terminals and 
forming a Cooper pair there. This is in contrast to the 
standard mechanism of (local) Andreev reflection"^ (AR) 
in which case two subgap electrons enter a superconduc- 
tor from the same normal electrode through the same 
interface. The phenomenon of crossed Andreev reflec- 
tion (CAR) manifests itself, e.g., in the dependence of 
the current II through the left NS interface of an NSN 
structure on the voltage Vr across the right NS interface. 
As a result, the non-local conductance Glr — dlL/dVji 
of an NSN device differs from zero and can be detected 
experimentally. Such experiments have recently been 
performed by several groupa^i^ providing a number of 
interesting observations some of which remain not fully 
understood. 

It is important to mention that CAR is not the only 
process which contributes to the non-local conductance 
Glr- Another relevant process is direct electron transfer 
(DET) between two normal terminals through the super- 
conductor. In the tunneling limit this process is nothing 
but the so-called elastic cotunncling (EC). It turned out^ 
that in the lowest order in tunneling the contributions 
from EC and CAR to Glr exactly cancel each other in 
the limit of low temperatures and voltages, i.e. the non- 
local conductance Glr should vanish in this limit. 

Note that this result^ is applicable only provided trans- 
missions of both NS interfaces remain small which is not 
always the case in the experiments. At higher transmis- 
sions processes to all orders should be taken into account 
and the contributions of DET and CAR do not anymore 
cancel each other. Hence, Glr does not vanish beyond 
the tunneling limit. In the case of ballistic electrodes 
a non-perturbative (in barrier transmissions) theory was 
recently developed by Kalenkov and one of the authors^iS,. 



This theory allowed to study the non-local conductance 
of NSN devices at arbitrary transmissions leading to a 
conclusion that CAR contribution to Glr vanishes in 
the limit of fully open NS barriers. This result might 
seem counterintuitive since ordinary (local) AR reaches 
its maximum at full barrier transmissions. In contrast, 
CAR is essentially a non-local effect which requires "mix- 
ing" of trajectories for electrons going between two nor- 
mal terminals with those for electrons going deep into a 
superconductor and describing the flow of Cooper pairs 
out of the contact area. Provided there exists no nor- 
mal electron reflection at both NS interfaces such mixing 
does not occur, CAR vanishes and the only remaining 
contribution to Glr in this case is one from DET. 

For completeness, let us point out that the exact can- 
cellation between EC and CAR contributions^ can also 
be violated by other means. One of them is simply to 
lift the spin degeneracy in the system. This can be 
achieved, e.g., by considering NSN structures with spin- 
active interfacesS or by using ferromagnets (F) as nor- 
mal metallic electrodeaiSiliii^. Experiments with FSF 
structures^ directly demonstrated the dependence of the 
non-local conductance Glr on the polarization of F- 
electrodes. 

Yet another way to avoid the cancellation between EC 
and CAR terms already in the tunneling limit is to in- 
clude interactions. This idea has been put forward in Rcf. 
[13I . The effect of electron-electron interactions on non- 
local conductance of NSN devices - in particular in the 
presence of disorder - is an interesting issue to be investi- 
gated further. Such investigation is, however, beyond the 
scope of the present paper. Here we only want to point 
out that interactions are not very likely to play the dom- 
inant role in the experiments^!^. This is because typi- 
cal resistances involved in these experiments were rather 
low and the corresponding dimensionless conductances 
strongly exceeded unity. Under such conditions the ef- 
fect of Coulomb interactions on AR is weaki^ii^ii^ and a 
similar situation with EC and CAR can be expected. In 
general, however, the combined effect of electron-electron 
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interactions and disorder on non-local properties of NSN 
devices can be important. 

In the absence of Coulomb interaction the effect of dis- 
order on non-local electron transport in NSN was recently 
considered in Refs. [millei. Brinkman and GoluboAsii 
employed the quasiclassical formalism of Usadel equa- 
tions and proceeded perturbatively in the interface trans- 
missions. They found that the proximity effect in the nor- 
mal electrodes in combination with disorder can strongly 
enhance both EC and CAR contributions to the non- 
local conductance. Duhot and Melin^^ argued that weak- 
localization-type of effects inside the superconductor may 
influence non-local electron transport in NSN structures. 
Morten et alJ^ considered a device with normal termi- 
nals attached to a superconductor via an additional nor- 
mal island (dot) and analyzed this structure within the 
framework of the circuit theory. In this paper we will 
extend and generalize this model by considering a super- 
conducting dot attached to one superconducting and two 
normal terminals as shown in Fig. 1. 

Our main goal is to study the combined effect of prox- 
imity and disorder inside the superconductor (dot). In 
addition, as it was demonstrated in experiments^"^, non- 
equilibrium effects, such as charge imbalance, inside a 
superconducting electrode may play a significant role. 
These effects will be included into our consideration too. 
We are going to show that the "peaked" temperature de- 
pendence of the non-local resistance Rlr observed in the 
experiments^!^ can be explained as a result of the compe- 
tition between charge imbalance and Andreev reflection. 
The crossover temperature between these two processes 
T* (defined in Eq. ([^5]) below) sets the position of the 
maximum in the dependence Rlr{T). 

In order to illustrate the main idea of our approach 
let us recall that the exact cancellation between EC and 
CAR terms^ occurs only at energies below the supercon- 
ducting gap while at higher energies (or in the normal 
state) CAR vanishes and EC remains the only relevant 
mechanism of electron transport. It is clear, therefore, 
that including the proximity effect due to the presence 
of normal electrodes immediately yields non-zero sub- 
gap density of states inside the supercoducting electrode 
which in turn should yield a decrease of CAR, thus elimi- 
nating its compensation by EC and leaving the non-local 
conductance Glr non-zero. This is precisely what we 
find. In order to correctly account for the above effects 
it is necessary to proceed non-perturbatively in tunnel- 
ing and consider interface conductances exceeding unity. 
Under these conditions the concept of elastic cotunnel- 
ing becomes irrelevant and it would be more appropriate 
to speak about direct electron transfer between two N- 
electrodes which includes processes of all orders in the 
interface transmissions. 

The structure of our paper is as follows. In Sec. 2 
we will introduce our model and outline the formalism 
to be used below. In Sec. 3 we will evaluate the Green- 
Keldysh functions of our system which will be used in 
Sec. 4 in order to derive the general expressions for the 
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FIG. 1: Superconducting quantum dot coupled to two normal 
and one superconducting leads. 

non-local currents in our NSN device. In the limit of low 
bias voltages these general results are further analyzed 
in details in Sec. 5 where we also illustrate their rela- 
tion to previous theoretical works and to experimental 
findings^'^. 

II. THE MODEL AND BASIC FORMALISM 

Below we will consider a chaotic superconducting 
quantum dot with the mean level spacing 5 connected 
to one superconducting (5) and two normal {L and R) 
massive electrodes by means of tunnel barriers. This 
structure is schematically shown in Fig. 1. The typ- 
ical size of the dot d is supposed to be sufficiently 
small, d ^; fo, where is the superconducting coher- 
ence length. Although we assume the channel transmis- 
sions of all three junctions are small, Tn'^ <^ 1 (here 
and below j — S, L, R), their dimensionless conductances 

gj — 2 TjP can take any (large) value provided the 
number of conducting channels is sufficiently large. The 
magnitudes of the superconducting order parameters in 
the dot and the electrode S are denoted respectively as 
A and Ag. The phase difference across the Josephson 
junction between the dot and the S-electrode is denoted 
by ip. 

Let us introduce the electron escape rate through the j- 
th junction T j — gj5 /A-k. Below we will demonstrate that 
transport properties of our system essentially depend on 
the parameters Fj/A. Our theoretical approach allows 
to obtain the exact solution of our problem applicable for 
all values of F^/A. 

The Hamiltonian of our system reads 

H — Hqd + Hl + Hr + Hs + Htl + Htr + Hts, (1) 
where 

Hqd^ J2 Cnt^Lcian + ^ (A4„d] _„ + C.C.) (2) 

n, a=T,i " 

is the Hamiltonian of an isolated quantum dot, 

Hl,R = ^ ekacl^.L,R^ka;L,R (3) 
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represent the Hamiltonians of the left and right normal 
leads, 



+ c.c. 



(4) 

is the Hamiltonian of the massive superconducting elec- 
trode and 



In our subsequent calculation we will make use of the 
fact that coupling of the n-th energy level of a chaotic 
dot to the leads does not depend on the number n, i.e. 
it remains the same for all levels. Hence, the effective 
level width Se = Tl + Tfi + Ts is also the same for all the 
dot levels. This observation enables us to first evaluate 
the 4x4 Green-Keldysh functions for each single energy 
level, then calculate its contribution to the current and 
afterwards perform a summation over all energy levels. 
This program will be accomplished below. 



define the tunneling Hamiltonians for the junctions j — 
S, i, R. 

Employing the Keldysh formalism we define non- 
equilibrium 4x4 Green-Kedysh function of rt-th energy 
level in the quantum dot: 



(7 F 




{dUtl)d\St2)) {T-^d^n{tl)d\jt2)) 



(6) 



III. GREEN-KELDYSH FUNCTIONS 

The Green-Keldysh functions of the dot Gin obey the 
Dyson equation 



where 



^-1 



E^in -A 

-E + in A 

-A E + £,n 

A Q ~E~£,„ 



, (9) 



G^ 



. / (rd]„(ii)'^'infe)> -{dUt2)d^iSh)) 

' \ {d\n{ti)dUt2)) {T-^d\jh)dat2)) 



and 



, tkm.itkn,j'^Gk,j{E)k. 



(10) 



nTd^r.{tl)d^n{t2)) -(di„(t2H„(tl)> 

V {d^n{ti)dln{t2)) {T-^d^n{tl)dln{t2)) ) ' 

/ {Td\^{h)d\^{h)) -{d\jh)d\jh)) 

V {d\jh)d]jt2)) (r-id]„(ti)4„(t2)) 



Note that here we do not introduce the off-diagonal ele- 
ments of the Green function Gnm ■ In the next section we 
will demonstrate that these off-diagonal elements vanish, 
Gnrn ~ for n ^ vfi^ provided the dot is fully chaotic 
and all channel transmissions are smaU, tI^^ < 1. The 
Green functions of the leads are defined analogously. 

The current through the left tunnel junction is ex- 
pressed as 



-GkAE)kGn{E)\. 



(7) 



Here A is the 4x4 diagonal matrix with the matrix ele- 
ments All — ~1, A22 = 1, A33 — 1 and A44 = —1. The 
functions Gn{E) and Gl{E) are the Fourier components 
of the Green-Keldysh functions for the dot and the left 
lead respectively. The currents across the right junction 
and across the Josephson junction between the dot and 
the superconducting electrode are defined analogously. 



is the self-energy of the j—th junction. 

In a chaotic quantum dot off-diagonal matrix elements 
of any operator between the m— th and the th energy 
levels tend to zero provided their energies are not too far 
from each other, |<^„ — ■Cml ^ D/d^ , where D is the diffu- 
sion coefficient for electrons inside the dot^. In addition, 
under these conditions the diagonal matrix elements do 
not depend on the level number n, i.e. (n| = const^. 
Hence, we obtain 



where 



ymn _ r y 



t,^\t,\'AY,GkAE)k- 



(11) 



(12) 



Then the Dyson equation for the Green-Keldysh function 
acquires the form 



-1 

n,qd 



G„ 



1. 



(13) 



Performing the summation over k in Eq. (|12p with the 
known expressions for the Green-Keldysh functions of the 
leads is straightforward. As a result, the self-energies for 
the junctions between the dot and the normal leads take 
the form 



R — 



^L,R ( a,Q{E±eVL,R) 



2i 









a,QiETeVL,R) 



(14) 
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where CTz is the PauU matrix and Q is the 2x2 matrix 
which reads 



^ (\-2n{E) 2n{E) 

" V 2 - ME) -1 + 2n{E) 



(15) 



Here n{E) — 1/(1 + e^/^) is the Fermi fmiction. 

The self-energy S5 for the Josephson jmiction between 
the dot and the S'-electrode, though somewhat more in- 
volved, is evaluated analogously. Combining the result- 
ing expression for Eg with Eq. (I14|) we obtain 



lip 



2 V -e 

.r(^) / 1 



TsFsjE) f E/As -e' 
{-e-''^ E/As 

2^ ) ^zQqd{E) 

Ts AsNsjE) f e'^ 
''2 E [ e""^ 



>a,QsiE), (16) 



where we denoted Fs(E) — ^^BjAs \e\) jntj-Qd^ced 
the density of states in the superconductor Ns{E) — 

\E\e{\E\-As) 



and defined 



TiE)^TL+TR + TsNs{E), 



(17) 



as the total escape rate of an electron from the dot 
through all three barriers. The Q-matrices in Eq. ([16]) 
read 



Qqd 



1 - 2nqd 2nqd 
2-2ngd -l + 2ngd J ' 



(18) 



where 



nUE) = ^fiE + eVL) + ^^n{E - cVr) 



2T{E) 



{Ns{E) + e{\E\ - As))n{E - eVs) 



+ ,^^{Ns{E)-e{\E\-As))n{E + eVs) (19) 

is the distribution function in the quantum dot, 

Q'gdiE, Vl,Vr, Vs) = Qgd[E, ~Vl.~Vr. -Vs) (20) 
and 

Qs(E) = {Q{E - eVs) + Q(E + eYs)\ /2. (21) 

We are now in a position to evaluate the Green- 
Keldysh function with the aid of the Dyson equation 
P^ . The derivation is facilitated by the normalization 
condition for the Q- matrices, = 1, as well as by the 
property 

i.Cj i 1 CJ i. ~\~ C/ j ■ 



After some algebra we finally arrive at the following ex- 
pression for the dot Green-Keldysh function: 

Gn{E,U = l{Gf,+G^)®^z 



+ i 



TsNsiE)As^Rf 



Gi 

2E " \ e 

r(g)^fl ^10^^^ 



^v ]Gn^QsiE)a, 



Q jGil® Qqd{E)d, 



(22) 



Here we defined the 2x2 matrix retarded and advanced 
Green functions 



qRA^^ ^ ip) - I GR,A{E,S.n,V) FR^A{E,£,n,^p) 



FR,AiE,^n,-^) GR,AiE,-^n,ip) 



where 

GrMe,^,^^) 

FR^A{E,t^) 



2As 



PrAe^Lv) 



, (23) 



A+^-^(Fs{E) 



Ns{E)As 
E 



PrAe,^,v) 



(24) 



are respectively the normal and anomalous retarded and 
advanced Green functions of the superconducting quan- 
tum dot and 

...(..,.,^(..I^..If)^ 

.2 , Tse'^ ,,,,^.Ns{E)As 



Fs{E)±i- 



E 



x(A.E^f^(f,W±.^«f^)). ,25, 

One can easily verify that the retarded and advanced 
Green functions are linked to each other by the standard 
relations Ga{E,£^,(p) = G*j^{E,£^,(p) and Fa{E,£,,(p) = 
F^{E,^-^). 

For completeness, we also present the self-consistency 
equation which controls the magnitude of the order pa- 
rameter in the dot: 



X [1 - n{E + eVs) - n{E - eVs)] 
X [cos2¥. GRiE,^„,^)G*j^{E,-^^,^) 

+ FR{E,^r.,^)F*RiE,^„.-v)] 
A 

+ 1 



I ^r(i?)Re[e^^G«(i?,C„,^) 
y.F*j,{E,^^,~^)\[l~2ngd{E)\. (26) 



where A is the BCS coupling constant. In general, the 
superconducting order parameter inside the dot should 
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be determined self-consistently with the aid of Eq. ((26)) . 
Here we avoid this complication and set A equal to a 
constant. This assumption is justified if, for instance, the 
coupling between the dot and the superconducting lead is 
much stronger than that between the dot and the normal 
leads, ^ Fl, Tr. If, in addition, we assume that both 
the dot and the superconducting lead are made of the 
same material, it would be appropriate to set A = A5 at 
all temperatures and sufficiently low bias voltages. 



IV. NON-LOCAL CURRENTS 

We now make use of the above general results and 
evaluate the currents across both NS interfaces of our 
device. Combining Eqs. ([7]) and ([12]) we express the 
current across the left interface in the form 

= IE / S*^' {[^l{E),A\G,.{E)). (27) 

An analogous formula is obtained for the current in the 
right jimction Iji. Substituting the results for the Green 
functions and self-energies derived in the previous section 
into the above expressions for the currents and setting 
Vs = 0, we obtain 



II = Ils{Vl) + ( 2^ + 1 ) Icar{Vl) + Idet{Vl) 



+ Idet{Vr) — Icar{Vr), 



(28) 



Ir - Irs{Vr) + ( 2^ + 1 ) Icar{Vr) + Idet{Vr) 



+ Idet{Vl)-Icar{Vl), 
where we defined 

/ls(^) = ^^E / dENsiE) 

n 

X 1\Gr{E,^^,^)\^ + \Fr{E,^^,^)\^' 



Re[G*R{E,^,,,^)FR{E,^„,-^)e'^] 



(29) 



2Ai 



X [n{E - eV) ~ n{E)] , 



(30) 



iRsiV) ^ f dENs{E) 

X \\GR{E,ir.,v)? + \FR{E,^^,^)\' 

Re[GR{E,Uv)F*R{E.£.n.^)e'- 



2Ag 

E 



X [n{E - eV) - n{E)] 



(31) 



Idet{V) = / dE\GR{E,i^) 



CAR. 



{V) 



X [n{E - eV) 
e^L^R 



TT 



(E)], (32) 
^ / dE\FRiE,^,,)\^ 

x[n{E~eV)-n{E)\. (33) 

Eqs. (l29 |) - ([33|) fully determine the currents across the 
left and the right NS interfaces and represent the central 
result of our paper. 

The current Idet{V) accounts for direct electron 
transfer between two normal terminals. This current dif- 
fers from zero also in the normal state of our system. 
In contrast, IcARiV) describes the contribution from 
crossed Andreev reflection which vanishes in the nor- 
mal limit. The contributions Ils and Irs contain terms 
which can be interpreted in a similar, though slightly 
more complicated manner since they originate from the 
Josephson junction between the superconductors and not 
from the NS interface. If, just for illustration, we put 
A5 = we immediately get Ils = {^s/^r){Idet + 
Icar) and Irs = (rs/TL){lDET + Icar)- 

We note that the possibility to decompose the currents 
Il,r, into the sum of partial currents each of which 
depending only on either Vl or Vr (but not on both) 
is due to the fact that the distribution function in the 
quantum dot nqd{E) (jl9p is represented as a linear com- 
bination of the distribution functions of the leads. This 
feature is similar to that of ballistic NSN devices^ii^. 



A. Normal state 

Let us analyze the above general expressions for the 
current. Considering first the trivial limit of a normal 
system A — A5 = we obtain 

TLn{E + eVl) + TRn{E - cVr) + Tsn{E) 

nqd = . (34) 

ti+ifi+is 

The current through the left junction is defined by a sim- 
ple formula 



1 



dE[nqd{E)-n{E + eVL)\. (35) 



The expression for Ir is similar. Evaluating the integral 
over E we obtain 



Vl 

RL ^RR I \ ^R 



where 



G 



N _ ^5 + Rr{L) 

LL{RR.) - j^^j^^ ^ ^^^^ ^ ^^^^ : 

Rs 



^LR — ^RL 



RlRs + RrRs + RlRr 



(36) 

(37) 
(38) 



In this limit both local and non-local differential conduc- 
tances remain voltage-independent. 
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In the experiments one often measures the non-local 
resistance 



Rlr 



dV, 



dlr 



G 



LR 



Ir=0 



GllGrr — G? 



LR 



From Eqs. ([36]) we obtain 



Rlr - R-S- 



(39) 



(40) 
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B. Charge imbalance 

Charge imbalance^ is a non-equilibrium phenomenon 
which is known to cause a number of interesting non-local 
effects in superconductors. This phenomenon is also of 
importance in connection with non-local electron trans- 
port in NSN hybrid structures discussed here. In par- 
ticular, it was argued^ that charge imbalance might be 
responsible for certain features of the non-local conduc- 
tance observed in experiments. Our approach allows to 
fully account for this phenomenon and its impact on non- 
local transport in the system under consideration. In this 
subsection we briefly illustrate the key physics associated 
with charge imbalance in our system. 

Just for the sake of illustration let us for a moment set 
A5 — and assume F^, Fj^, F5 ^ A. In this regime the 
current II ([^5]) across the left junction takes the form^^ 



II = 



eRL 



dE '^'^i^':^) \n{E) n{E + cVl)] 



cRl 



y/E^ - A2 
dE9i\E\~A) fiqd{E), 



(41) 



where hqd{E) is the asymmetric part of the distribution 
function responsible for charge imbalanceS^. For the sys- 
tem under consideration nqd{E) reads^^ 



lqd{E) 



_ A2 nq^(^E) + nqd{-E) - 1 



\E\ 



(42) 



and nqd{E) is given by Eq. (|34|) . Then for local and 
non-local zero bias conductances one obtains 



G 



1 



LL{RR) 



dE 



0{\E\-A) 



\E\ 



Rl{r)J '~ 4Tcosh2J. \VE^-A^^ 
RRiL)Rs VE^-AA 



RlRr + RlRs + RrRs 



\E\ 



Glr. — G 



N 
LR 



dE- 



\E\-A) y/W~^ 



4rcosh 



2 _E_ 

2T 



\E\ 



,(43) 



(44) 



Note that in this regime the non-local conductance is 
solely due to charge imbalance being fully determined by 
the second term in Eq. (|4ip . At low temperatures we 
have GiL, Gfljj, Gifl oc e~^l'^ . Hence, in the situation 
considered in this subsection at T the non-local 
resistance ([Ml) should diverge as Rlr oc i?5e^/^. 



FIG. 2: Normalized local Gll{Vr)IGIj^ and non-local 
GhRiyR.) /GIji differential conductances as a function of ap- 
plied voltage Vr at low temperature T ^ A and at different 
tunneling rates Tl = Fn. Here we set Ag = A and Fg = lOA. 
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FIG. 3: The same as in Fig. [2] for different values of the dot 
gap A. Here we set Tl=Tr = 0.02As and Vs = O.lAs. 



C. General case 

Now let us return to the case A5 ^ 0. With the aid 
of Eqs. ([5^ . ([55)1 we determine the differential conduc- 
tance Glr which is presented in Fig. [2] as a function 
of applied voltage Vr for A = A5 and different values 
of the tunneling rates F^ = F^ as compared to A. We 
observe that at subgap voltages cVr < A the magni- 
tude of the normalized non-local conductance Glr/Gn 
increases with increasing {Tl + T]i)/A. Such dependence 
is quite natural because exactly the same ratio controls 
the strength of the proximity effect in our system. As we 
have already discussed, with increasing value of the ra- 
tio (Fi -|- Tfi)/ A the proximity-induced subgap electron 
density of states increases, the difference between DET 
and CAR contributions grows and, hence, Glr becomes 
bigger. 

Fig. [3] illustrates the effect of the dot order param- 
eter A on the non-local conductance. At high volt- 
ages eVfl, > A5 we recover the normal state value ([55)1 
while at intermediate values A < eVfl ^ A5 we find 
Glr ~ ^/{Rl + Rr)- For eVn < A the conductance 
Glr progressively increases with decreasing ratio A/ As 
and eventually reaches the maximum in the limit A = 
in which case the resultsi^ are reproduced. 

In order to demonstrate an important difference in the 
low voltage behavior of Glr for superconducting and 
normal quantum dots in Fig. [3] we deliberately chose 
small values of tunneling rates Tl,r ^ A. We observe 
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that in the superconducting case Glr essentially van- 
ishes at eVfi < A in accordance with Ref. 0, while for 
normal quantum dots^^ Glr remains non-zero even at 
Vb, ^ 0. 



Considering a big metallic quantum dot we can replace 
the sum over energy states by the integral ^„ ~^ \ ! 
Making use of the relation between the tunneling rates 
and the the junction resistances, Tl^h — 5 /2e^RL^ji, we 
reproduce the standard result 



V. ZERO-BIAS CONDUCTANCES 

Let us now consider the behavior of the conductance 
matrix in the limit of low voltages in more details. In the 
zero bias regime currents flowing through the system are 
low, i.e. we can set if — Q. Similarly to Ref. 7 at low 
voltages the expression for current 1^ can be split into 
three different contributions 

II - GaVl + Gdet{Vl + Vii) + Gcar{Vl ~ Vb). (45) 



G 



DET 



Rl + Rr. ' 



(51) 



i.e. in this case DET contribution simply reduces to the 
Ohm's law. 

Next we put T b. ~ and consider a superconducting 
dot coupled to one normal and one superconducting lead. 
In this case one trivially gets Gdet = Gcar = 0. Pro- 
vided ri,rs ^ A the dot can be viewed as a point-like 
scatterer with the following set of transmission probabil- 
ities (cf. Eq. ([501)1 



Here Ga is (local) Andreev conductance of the left NS 
barrier, Gdet and Gcar are respectively DET and CAR 
contributions to the zero bias conductance matrix. 



A. Zero temperature limit 

In the limit of zero temperature T ^ from Eqs. ([29]) - 
([55)1 we obtain 



2Tl (a2 + TsA + if 



G 



DET 



e 

TT 



E 



(rL+rB)^+r| 



TsA 



(46) 



-2 , {rL+rny 

4 



,2 , (rL+Tny+rl 

^ ^ 4 



r.,A 



r(47) 



TrT 



S 



£2 



(52) 



We note that, although the channel transmissions of NS 
interfaces remain small, effective transmissions f„ are not 
necessarily small. The Andreev conductance in this limit 
becomes 



Ga 



2tt ^ 



5? 



(53) 



One can verify that this expression can be cast to the 
familiar form^iiS^ 



(54) 



In a general case of metallic quantum dots one can 
perform the summation over £„ in Eqs. (|46ll48p and arrive 
at the following explicit expressions 



Gcar 



L4 R, 



A2 -t- 



EcA 



..(48) 



Let us analyze the above expressions in different phys- 
ical limits. We first put Eg = and A = 0, i.e. we 
consider a normal quantum dot isolated from the su- 
perconducting electrode. Then we obviously find G a = 
Gcar. — 0, while for Gdet we obtain 



Ga 



e2 g+(.9L+gfl)V2 

^DET = ^9LgR -j^^ , 



B 

GcAR = -^gLgRj^- 



(55) 



(56) 



(57) 



G 



DET 



TrT 



Li R 



(49) 



Comparing this expression to the Landauer formula we 
immediately conclude that each energy level of the dot 
effectively corresponds to one conducting channel with 
transmission 



ErE 



Lr B 



(50) 



where 



and 



B 



St- 



as- 



is, 

4 



<52 4 



47rA 



-gs- 



(58) 



(59) 



In the limit A ^ our results for Gdet and Gcar 
reduce to the corresponding expressions derived in Ref. 
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llSl for the normal quantum dot. At the same time, our 
result ([55]) for the Andreev conductance Ga (for A 
0) turns out to be 4 times bigger than the analogous 
expression^^. This difference is supposed to be due to a 
different definition of the Andreev conductance employed 
in Ref. HI. 

Combining the above results for Gdet and Gcar we 
immediately arrive at the zero temperature linear non- 
local conductance Glr = Gdet — Gcar for our device. 
It reads 



G 



LR 



9L9R{gL+ qrY 



IGtt 



16Tr2A2 



{gh+gRY+gl 



r, 47rA 

.95 — 



3/2 



(60) 



This expression demonstrates again why the lowest order 
perturbation theory in barrier transmissions- yields zero 
non-local conductance at T = 0. This perturbation the- 
ory applies in the weak tunneling limit gL,R ^ 1- The 
result (|60p . however, contains only higher order terms in 
barrier transmissions whereas the contribution oc gL9R. 
should vanish. This situation is qualitatively similar to 
that of NSN structures with ballistic electrodes^i^. We 
would also like to emphasize that the exact cancellation 
of Gdet and Gcar in the lowest order in ghdR holds for 
any gs and does not require taking the limit gs ^ oo. 
This is in contrast to the case of normal quantum dotsi^ 
in which Glr was found to vanish only for gs oo. 

At small tunneling rates rL,r_R,rs < A Eq. ((SD|) 
reduces to 



Glr — 



(167r) 



9LgR{9L + gR)'^S^ 
- A3 



(61) 



In the limit of a bulk metal 6^0 (though d ^ ^o) the 
proximity effect becomes unimportant and the non-local 
conductance Glr (pT|) vanishes already to all orders in 
gL,R- 

Finally, we present the exact expression for the zero 
temperature non-local resistance Rlr- It reads 



R 



LR 



2r. 



(( 



A- 



4 



3/2 



Rs 



(2(A + ^#) 



4 



-{TL + TRy 



.(62) 



In the limit Ts ^ AjTlj^r we get Rlr = Rs, i-e. 
in this limit the non-local resistance just coincides with 
its normal state value. For A ^ ri,r/i,rs we obtain 
Rlr = 5/4e^A « Rg. 



B. Non-zero temperatures 

Finally let us briefly discuss the effect of temperature 
on zero bias conductances of our system. Combining Eqs. 
- ([55]) and (gSl) we obtain 



G 



LL 



dE 



\E\ei\E\-As) 



10" 
10-' 
10-^ 



+ 

£ IxlO"' 
Ixio" 
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(b) 
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FIG. 4: (a) temperature dependence of local Gll (|63[) and 
non-local Glr (|64I) zero- bias conductances. The parameters 
were chosen as follows: A = As, Tl = = O.IA(O), 
Fs = lOA(O); (b) non-local zero bias resistance Rlr (p9|) 
normalized by its normal state value (for the same parame- 
ters). 



\GRiE,Uf + \Fn{E,U\' 



4Tcosh^{E/2T) 

Re[GR{E,^r.)F*n{E,Cn)] 



2As 

E 



Li R 



E 



dE 



ATcos\i\E/2T) 



X { \Gr{E,U)\^ + ( E£ + 1 ) \FR{E,i^)\^ \ ,(63) 



Glr. — 



e'TLVR 



dE 



\GR{E,i,,)\''~\FR{E,i,,)\ 
4Tcoah^{E/2T) 



(64) 



Substituting the expressions for the Green functions ([23])- 
(1^5)1 into the above equations we arrive at the final results 
for zero-bias conductances at non-zero T. These results 
are illustrated in Figs. 4 and 5. 

Fig. 4 shows the temperature dependence of both lo- 
cal and non-local zero bias conductances ((63|) , ((64|) along 
with that for the non-local resistance ([39| . The conduc- 
tances Gll and Glr. decrease monotonously with de- 
creasing temperature. The temperature dependence of 
the non-local resistance Rlr is, on the contrary, non- 
monotonous. At temperatures just below Tc the re- 
sistance Rlr first slightly decreases but then it starts 
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0,0 0,2 0,4 0,6 0,8 1,0 1,2 



FIG. 5: Temperature dependence of zero-bias non-local con- 
ductance Glr for different values of Fg. Here we set A = As, 
Fl = Fjj = 0.3A(0). 

growing exponentially with decreasing T due to charge 
imbalance effects, as it was explained in Sec. 3b. Such a 
tendency persists down to the crossover temperature 

T*-7-^+rZ^ ri,rfl«2A + rs (65) 

at which the non-local resistance reaches its maximum. 
Below T* , AR contribution starts dominating over that 
caused by charge imbalance. For this reason at T < T* 
Rlr drops sharply and then at T ~ T* /2 saturates to 
its zero temperature value ([5^ . as it is seen in Fig. 4. 

Note that qualitatively the same behavior of the non- 
local resistance was recently observed in experiments^!^, 
cf., eg.. Fig. 3 in Ref. H. Though a detailed quantitative 
comparison between our theoretical predictions and the 
experimental results^i^ is rather difficult to perform due 
to different geometry of the model employed here, we 
believe that our theory correctly describes the physical 
origin of the peak in the temperature dependence of the 
non-local resistance observed in Refs. 011. 

For completeness, in Fig. 5 we present the dependence 
Glr{T) at different values of Fg. As temperature de- 
creases below the critical temperature Tc the conduc- 
tance Glr{T) drops sharply below its normal state value 
(155)) and at T <C Tc it saturates to the zero-temperature 
value which essentially depends of the relation between 
^L,R,s and A. We observe that for given Tl.r this 
value decreases with increasing coupling Fg between the 
dot and the superconducting lead. This tendency is ex- 
plained by the fact that CAR becomes progressively more 
pronounced with increasing tunneling rate Fg. 



which consist of a superconducting chaotic quantum dot 
(with typical size d ^ ^q) attached to one superconduct- 
ing and two normal reservoirs, as it is shown in Fig. 1. 
By varying the tunneling rates between the dot and the 
electrodes Tl^r^s (which play the role of effective Thou- 
less energies for electrons in the part of a superconduct- 
ing electrode directly attached to normal leads) one can 
cover a number of different physical situations and limits 
and illustrate the relation to the models considered by 
other authors. 



Our analysis is employed within the general Keldysh 
formalism which fully accounts for non-equilibrium ef- 
fects and disorder in the superconducting terminal (dot). 
Our theory allows to go beyond perturbation theory in 
dimensionless conductances between S- and N-electrodes 
gL,R and derive a general expression for the conductance 
matrix which remains valid in both weak and strong tun- 
neling limits. This result enables one to study and com- 
pare relative contributions to the non-local conductance 
provided by the competing processes of direct electron 
transfer (DET) and crossed Andreev reflection (CAR). 
We demonstrated that at low energies these contribu- 
tions do not cancel each other beyond the weak tunnel- 
ing limit. This is the result of the proximity effect: Cou- 
pling to normal electrodes yields non-zero subgap den- 
sity of states inside the superconducting dot which in 
turn causes a decrease of the CAR contribution to the 
non-local conductance Glr- On the contrary, increas- 
ing coupling between the dot and the superconducting 
electrode increases CAR and, hence, decreases Glr- 



Our theory allows to investigate the effect of charge 
imbalance on non-local electron transport in NSN de- 
vices. We argued that temperature dependence of the 
non-local resistance Rlr of such devices is determined 
by the competition between charge imbalance and An- 
dreev reflection. The contribution of the former process 
dominates over that of the latter at T > T* (where T* is 
the crossover temperature defined in Eq. (|55|) ) causing 
an increase Rlr{T) with decreasing T. In contrast, at 
lower temperatures AR dominates and Rlr{T) decreases 
as T becomes lower. As a result, the dependence Rlr{T) 
acquires a pronounced peak at T ~ T*, see Fig. 4. This 
behavior was observed in recent experiments^i^. 



VI. CONCLUSIONS 

In this paper we developed a microscopic theory of non- 
local electron transport in three-terminal NSN structures 
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